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ABSTRACT 

Fermi-LAT has recently detected two gamma ray bubbles disposed symmetrically 
with respect to the Galactic plane. The bubbles have been suggested to be in a quasi- 
steady state, inflated by ongoing star formation over the age of the Galaxy. Here we 
propose an alternative picture where the bubbles are the remnants of a large-scale 
wide-angle outflow from Sgr A*, the SMBH of our Galaxy. Such an outflow would 
be a natural consequence of a short but bright accretion event on to Sgr A* if it 
happened concurrently with the well known star formation event in the inner 0.5 
pc of the Milky Way ~ 6 Myr ago. We find that the hypothesised near-spherical 
outflow is focussed into a pair of symmetrical lobes by the greater gas pressure along 
the Galactic plane. The outflow shocks against the interstellar gas in the Galaxy 
bulge. Gamma-ray emission could be powered by cosmic rays created by either Sgr A* 
directly or accelerated in the shocks with the external medium. The Galaxy disc 
remains unaffected, agreeing with recent observational evidence that supermassive 
black holes do not correlate with galaxy disc properties. We estimate that an accreted 
mass ~ 2 x 1O 3 M0 is needed for the accretion event to power the observed Fermi- 
LAT lobes. Within a factor of a few this agrees with the mass of the young stars 
born during the star formation event. This estimate suggests that roughly 50% of the 
gas was turned into stars, while the rest accreted onto Sgr A*. One interpretation 
of this is a reduced star formation efficiency inside the Sgr A* accretion disc due to 
stellar feedback, and the other a peculiar mass deposition geometry that resulted in a 
significant amount of gas falling directly inside the inner ~ 0.03 pc of the Galaxy. 
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1 INTRODUCTION 

1.1 Sgr A* - the SMBH of the Milky Way 

Sgr A* is the supermassive black hole (SMBH) in the nucleus 
of our Galaxy. Its mass M bh ~ 4 x 10 6 M Q (Schodel et al. 
2002, Ghez et al. 2005, 2008) makes it directly comparable 
with SMBH in other galaxies. The Soltan relation (Soltan 
1982) implies that most of the mass of these black holes was 
gained through luminous accretion. Yet by comparison with 
active galactic nuclei (AGN) Sgr A* is famously dim. It is 
spectacularly faint both in X-rays (less than ~ 10~ Z/Edd, 
where LEdd ~ a few x 10 44 erg s _1 is its Eddington luminos- 
ity; Baganoff et al. 2003) and in the near infrared (Genzel et 
al. 2003), prompting suggestions of a radiatively inefficient 
accretion flow (Narayan 2002, and references therein). Cur- 
rently, Sgr A* appears to be fed by accretion of gas captured 
(Cuadra et al. 2006) from the winds of the young massive 
stars populating the inner ~ 0.5 pc of the Galaxy (Paumard 



et al. 2006). However, X-ray reflection nebulae suggest that 
Sgr A* might have been much brighter a few hundred years 
ago, with luminosity of a few xlO 39 erg s _1 (e.g., Revnivt- 
sev et al 2004, Ponti et al 2010). This may reflect variations 
in the wind feeding rate of Sgr A* caused by changes in the 
stellar orbits of the most important wind-producing mas- 
sive stars (Cuadra et al. 2008), or longer time scale feeding 
events from a few pc-scale molecular gas reservoirs (Morris 
et al. 1999). 

Sgr A* is also famous as the site of a recent (~6x 10 6 
yr ago) star formation event in one and perhaps two stellar 
discs (Levin & Beloborodov 2003, Genzel et al. 2003, Lu et 
al. 2009) on scales of ~ 0.03 - 0.5 pc from the SMBH. The 
observed (e.g. Paumard et al. 2006) and theoretically con- 
strained (Nayakshin et al., 2006) mass of the young stars is 
around a few times 1O 3 M0, perhaps even 1O 4 M . Signifi- 
cantly, there is currently no trace of even a remnant gaseous 
disc near Sgr A* (Cuadra et al. 2003, Paumard et al. 2004). 
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This led Nayakshin & Cuadra (2005) to question whether 
Sgr A* failed to become a quasar because this recent star 
formation event consumed nearly all the available gas in the 
central parsec of the Milky Way. They noted that this could 
be constrained with future observations: "a past bright AGN 
phase should also leave a hot buoyant radio bubble in the 
Milky Way halo". 

1.2 The Fermi— L AT gamma— ray lobes 

The recent Fermi-LAT observations by Su et al. (2010) show 
that the Milky Way has a pair of gamma-ray lobes, sym- 
metrical about its dynamical centre (Sgr A*) and about the 
Galactic plane. The lobes extend ~ 5 kpc from the plane, 
but have a narrow (d ~ 100 pc) waist along it. The limbs 
of the lobes coincide with the extended structure seen in 
medium-energy X-rays by ROSAT (Snowden et al. 1997). 
The lobes have gamma-ray luminosity L 7 ~ 4x 10 37 ergs -1 , 
and their total energy content is at least ~ ^q 54 - 55 er g 

Su et al. (2010) considered numerous physical processes 
that could give rise to the bubble structure and provided 
a constraint that if they are older than a few x 10 6 yr, the 
gamma-ray emission must be powered by ions rather than 
electrons due to a short cooling time of the latter. Crocker et 
al. (2011) and Crocker & Aharonian (2011) detailed these ar- 
guments further and suggested that the emission is powered 
by Cosmic Ray (CR) protons rather than electrons. They 
further consider a quasi-steady state model in which the 
CR protons are continuously injected by supernova explo- 
sions. CR protons and heavier ions are then trapped inside 
the bubbles for approximately the age of the Galaxy. 

Alternatively, the Fermi-LAT lobes could be a more re- 
cent feature. For example, Mertsch & Sarkar (2011) argue 
that the spectral and morphological details of the emission 
are incompatible with hadronic radiation, leaving electrons 
as the main energy source. In that case, the mechanism of 
inflating the bubbles is then unlikely to be of star forma- 
tion origin. One would require ~ 10 5 recent Type II super- 
novae to provide the energy content of the bubbles, which 
is far higher than can be realistically expected from the in- 
ner ~ 100 pc. Cheng et al. (2011) thus argued that the 
bubbles are inflated by episodic Sgr A* activity caused by 
tidal disruptions of stars passing too close to Sgr A* . Guo &: 
Mathews (2011) performed hydrodynamical numerical sim- 
ulations of jets launched by Sgr A* and showed that the 
Fermi-LAT observations are qualitatively consistent with 
their simulations if the jets were launched ~ 1 — 2 Myrs ago. 

1.3 Sgr A* feedback: when and how? 

In this paper we shall argue that Sgr A* is a very natu- 
ral candidate for the source of the energy that inflated the 
gamma-ray lobes. As noted above, the Galactic Centre un- 
derwent a peculiar star formation event localised to the inner 
0.03 - 0.5 pc about 6 Myr ago (Paumard et al. 2006). Thus, 
a plausible scenario is that not all of the gas deposited into 
the central pc of the Milky Way went into making the young 
stars, and a fraction of it was accreted by Sgr A*, as found 
in realistic simulations of the process (e.g., Bonnell & Rice 
2008, Hobbs & Nayakshin 2009). Thus Sgr A* is likely to 
have had a short but very bright quasar phase concurrent 
with the star formation event ~ 6 million years ago. 



We further argue that the observed highly symmetrical 
lobes are unlikely to have originated from a jet outflow. To 
obtain the qualitative agreement with the observed shape 
of the lobes, Guo & Mathews (2011) directed their jets per- 
pendicular to the plane of the Galaxy. We believe this would 
be unlikely. Radio surveys show that jet directions are com- 
pletely uncorrelated with the large-scale structure of the 
host galaxies (Kinney, 2000; Nagar & Wilson, 1999). Fur- 
themore, the observed orientations of the stellar discs in the 
central pc of the Galaxy (see Paumard et al. 2006) are in- 
clined at very large angles to the Galactic plane. The jets 
are likely to be fed by gas discs oriented similarly to the 
stellar discs. We would therefore expect that accretion of 
gas onto Sgr A* ~ 6 million years ago would result in jets 
directed at very large angles to the symmetry axis of the 
lobes, contradicting observations. 

In contrast, a symmetrical pair of lobes with a nar- 
row waist along the galaxy plane is natural if an isotropic 
outflow from near the black hole encounters higher gas 
densities along this plane than perpendicular to it. Near- 
spherical outflows like this are a direct consequence of super- 
Eddington disc accretion (Shakura & Sunyaev 1973; King & 
Pounds 2003) and offer a plausible explanation for the M—a 
relation (Silk & Rees 1998; King 2003, 2005). 

The paper is structured as follows. We first discuss the 
simpler and better understood quasi-spherical AGN outflows 
in §2, and then we consider the more complicated case of 
the present day Milky Way nucleus in §3. The implications 
of the quasar outburst for the observed gamma-ray lobes 
are elucidated in §4, while §5 spells out ramifications for 
the poorly understood problem of star formation versus gas 
accretion in the central parsecs of AGN. We note that our 
approach here is to try to reproduce the main energetics 
of the lobes and their morphology rather than to produce 
detailed spectral models. We assume that Sgr A* outflow 
either carries with it CR protons created near the black hole, 
or that the CR protons are accelerated on shock fronts where 
the outflow runs into the interstellar medium. 



2 SPHERICAL OUTFLOWS 

In regions close to the black hole, the AGN outflows are re- 
vealed through blueshifted absorption lines in X-ray emis- 
sion (Pounds et al. 2003a, b; King 2010a). Tombesi et al. 
(2010a, b) show that they are present in more than 35 per- 
cent of a sample of over 50 local AGN, and deduce that 
their solid angles are large (certainly > 0.6 x 2n, and prob- 
ably greater). The observed absorption columns imply that 
in many cases the outflows are quite recent (few years) , sug- 
gesting that outflows are an almost ubiquitous feature of 
central black hole activity (King 2010b). 

Although supermassive black holes (SMBHs) in galaxy 
centres frequently accrete at the Eddington rate, accretion 
at significantly higher rates requires extreme conditions (cf 
King, 2010a). Accordingly we consider cases where the ac- 
cretion rate far from the SMBH only mildly exceeds A^Edd, 
and both the central accretion rate and the outflow rate 
M out are ~ M E dd- Then the outflow has scattering optical 
depth ~ 1, and the photons driving it typically scatter only 
once before escaping. The front-back symmetry of electron 
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scattering means that the outflow momentum must be of 
the same order as the original photon momentum, i.e. 

^Edd 



M out v 



(1) 



so that the outflow velocity v ~ r/c, where 77 ~ 0.1 is the 
accretion efficiency (e.g. King & Pounds 2003; King 2010a). 
The wind flows with essentially constant velocity v until it 
shocks against the interstellar gas of the host galaxy, driv- 
ing a second shock outwards into this ambient medium and 
sweeping its up into a shell. A simple representation of the 
interstellar density is the isothermal distribution 



p(R) 



f a o 2 



(2) 
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where f g is the gas fraction, and a is the velocity dispersion. 
The average cosmological value of f g is f c — 0.16. 

Within this model, then, in galaxies with large a ^ 150 
km s _1 , Eddington outflows tend to sweep the vicinity of the 
hole clear of gas of density ([2]) and prevent further accretion 
and growth, establishing the M— a relation for the black hole 
mass (King 2003; 2005). At smaller values of a, any effect of 
this kind is outdone by the effects of mass loss from nuclear 
star clusters. These sweep out the gas (McLaughlin et al., 
2006; Nayakshin et al., 2009) and establish an offset M- a 
relation between the total cluster mass and the bulge veloc- 
ity dispersion (Ferrarese et al. 2006, Fig. 2, middle panel). 
The Milky Way is probably a member of this star-cluster 
dominated class of galaxies, and indeed its SMBH mass lies 
significantly below the value predicted from the M-a rela- 
tion (see, e.g. Greene et al. 2010, Fig. 9). 

The double shock pattern caused by the impact of an 
Eddington outflow on the host interstellar medium must 
move radially outwards from the vicinity of the black hole. 
The nature of this motion depends crucially on whether or 
not the shocked wind cools within the flow time. If cooling 
is effective, most of the energy injection rate 



E = ~ M out v 2 



2 2 

— —Mout — — -t^Edd 



(3) 



is lost to radiation, and only the ram pressure of the outflow 
is communicated to the host ISM. This is a momentum- 
driven flow. If instead the flow does not cool, the shocked 
wind expands adiabatically, doing PdV work against the 
swept -up interstellar medium. This is an energy-driven flow, 
which expands faster through the ISM than a momentum- 
driven one. 



3 FERMI— LAT LOBES AS QUASI-SPHERICAL 
OUTFLOWS 

The gamma-ray lobes observed by the Fermi-LAT instru- 
ment are very wide features that we shall first consider ap- 
proximately quasi-spherical. For the present day Milky Way 
and directions well out of the Galactic plane, we expect f g 
to be significantly less than f c , so we parametrize f g as 
f g = 1.6 • 10~ 3 /o.oi, where /0.01 ~ 1 is a dimensionless free 
parameter of the model. 

We now ask if the shocked gas cools in conditions appro- 
priate for the outburst. The outflow speed v ~ 0.1c implies a 
shock temperature T s = (3m p /16k)r) 2 c 2 ~ 1.6 x 10 10 K. This 
is much higher than the Compton temperature ~ 10 7 K of 



the SMBH accretion flow, so when the shock is sufficiently 
close to the hole, Compton cooling by the central radia- 
tion field is very effective and enforces momentum-driven 
flow. As the shock reaches a critical radius -R en the radi- 
ation field becomes too dilute to cool it. Also, the shocked 
wind has far too low a density to cool effectively by atomic or 
free-free processes, so the flow becomes energy-driven (King 
2003; King et al. 2011). For the parameters of Sgr A* (mass 
Mbh — 4 x 10 6 M Q , velocity dispersion a ~ 100 kms' 1 ) the 
transition to energy-driven flow occurs at a radius 

R cn ~ pc (4) 



(cf eqs 8 - 10 of King 2003). 

Even at the cosmological gas fraction (/0.01 = 100) 
the estimate is so small compared with the size of the 
gamma-ray lobes that we can regard the outflow as es- 
sentially always energy-driven in directions away from the 
Galactic plane. In an energy-driven outflow, the shocked 
wind density driving the expansion is always much lower 
than the density of the swept-up interstellar medium out- 
side it. This makes the shock interface inherently Rayleigh- 
Taylor unstable (cf King 2010b). The hot shocked gas mixes 
with cool dense interstellar gas throughout the flow in direc- 
tions away from the Galactic plane. This mixture is clearly a 
promising site for gamma-ray emission. Within the Galac- 
tic plane the gas density is far higher, and we expect lit- 
tle expansion (see also §4). This kind of outflow thus natu- 
rally produces the main qualitative features of the Fermi- 
LAT gamma-ray map: extensive gamma-ray emitting lobes 
placed symmetrically on each side of the Galactic plane, with 
a narrow waist in the plane. 

Energy-driven outflows rapidly attain a constant speed 
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(5) 



in the bulge of a galaxy (King 2005) . Here the factor b < 1 
allows for some collimation of the outflow, M a is the pre- 
dicted value of the SMBH mass in the Milky Way from the 
M — a relation and M ~ 0.2M CT is the mass of Sgr A*. 

For the rest of this Section we model the outflow away 
from the disc plane as a sector of a spherical flow. If the 
Eddington accretion phase lasts for a time t acc , the shock 
reaches radius 

R ~ v e t acc (6) 

when the quasar phase ends. However the shocked wind gas 
is able to drive further expansion, which finally stalls at a 
radius 



Rstall — -Ro — — tacc 

a a 
after a time 



i .... ,ii u.;> [ — ) / . 
(King et al. 2011) so that 

-Rstall 



tstall 



2a 



(7) 



(8) 



(9) 



We now apply these considerations to our Galaxy, and 
in particular the gamma-ray lobes. If the outflow producing 
the observed lobes had stalled, we would have -Rstaii ~ 5 kpc, 
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which from ((9| requires t s tall — 15 Myr. This would mean 
that the outflow was produced well before the last accretion 
event in the Galactic Centre, which appears unlikely. 

If instead we assume that the gamma-ray lobes were 
produced in this event, we must conclude that the energy- 
driven outflow is still proceeding, with a mean velocity (v) — 
1000 kms -1 over its lifetime. This is lower than the shell 
velocity during the quasar phase, which might have been as 
large as v e — 1600 km s _1 (cf. the figures in King et al. 2011). 
This is compatible with ([5]) if /o.oi ~ 1, i.e. if the gas fraction 
in the lobes is about 1% of the cosmological value. Requiring 
tstaU > 6 Myr in <[Sj) now gives £ acc >5 x 10 4 yr. At the 
Eddington accretion rate ~ 4 x 10~ 2 M0yr _1 appropriate 
for the mass of Sgr A* this gives the total mass accreted by 
the black hole during the quasar phase as AM > 2 x 10 3 Mq. 
This is comparable to the total expected if the hole accreted 
the disc mass 

M disc ~ |m bh ^ 8000 Mq (10) 

within the self-gravity radius where the ring of young stars 
formed (cf eqs 7, 12 of King & Pringle, 2007). This estimate 
is also consistent with the results of Nayakshin & Cuadra 
(2005). 

In this picture the mass of wind expelled from the vicin- 
ity of the hole must be ~ AM. Almost all of its kinetic 
energy is retained by the outflow. This energy is of order 

2 

E lobcs ~ ?L AMc 2 > 2 x 10 55 erg, (11) 

somewhat above the minimum required by observation. At 
the current gamma-ray luminosity the lifetime is ~ 10 10 yr, 
but there may be other losses of course. We conclude that 
the properties of the gamma-ray lobes are consistent with 
their production in a short phase of Eddington accretion 
about 6 Myr ago. 



which shows convincingly that the outflow from Sgr A*, even 
in its full "quasar" mode, is not strong enough to disperse 
the CMZ since the latter is simply too massive. This con- 
clusion is reinforced by the fact that there is also atomic 
and ionised gas in the region of the CMZ disc that would 
increase W C ms> further. 

Another way to come to the same conclusion is through 
estimation of the gas density in the midplane of the CMZ, for 
which we infer p C mz ~ 5 x 10~ 22 g cm -3 with the parameters 
mentioned above, whereas the density of gas which could be 
driven away by a SMBH outflow, for a SMBH obeying the 
M — a relation, is given by equation [2] and is ~ 10 -22 g 
cm" 3 at R = 200 pc and a = 100 km s" 1 . 

We therefore conclude that the outflow along equatorial 
directions stalls. As argued above, the outflow should then 
thermalise and expand away from the plane of the symmetry, 
i.e., the Galactic plane. This would naturally explain the 
two-lobe structure of the Fermi-LAT bubble. 

We also note that a relatively geometrically thin dis- 
tribution of the molecular gas along the Galaxy plane in 
the CMZ justifies our assumption that gas mass and density 
in the direction significantly away from the plane is low. In- 
deed, the density at height z ~ R away from the midplane is 
~ exp[~(R/H) 2 ] that in the midplane. Even at H/R = 0.3 
this factor is about 10 -3 . We therefore estimate that the 
outflow should have an opening angle larger than 45°, and 
realistically in the range of ~ 60 — 70°. 

Thus although outflows ultimately control black-hole 
growth, and materially affect the Galaxy bulge, as shown 
by the M — a relation, they cannot disperse the Galaxy 
disc. This fits very naturally with the recent conclusion by 
Kormendy, Bender & Cornell (2011) that SMBHs do not 
correlate observationally with host galaxy discs. 



4 THE NARROW WAIST OF THE 
GAMMA-RAY LOBES 

So far we considered the outflow in a spherical geometry. 
However, it is well known that the central ~ 200 pc of the 
Galaxy host as much as ~ 10% of all molecular gas of the 
Galaxy in a flattened, presumably disc-like, configuration 
(Morris & Serabyn 1996). We shall now show that this fea- 
ture, called the Central Molecular Zone (CMZ), could not 
have been significantly affected by the hypothesised outflow 
from Sgr A*. 

The mass of the molecular gas in the zone is M cmz ~ 

5 x 1O 7 M . Its weight is 

GM C 

nc Ukmz )- 1, ' J cmz 

2M cmz cr 2 , . 

Wcmz ~ ^ '- = — , (12) 

where M enc (_R cmz ) is the mass enclosed within radius i? cmz ~ 
200 pc. The outward force (momentum flux of the outflow 
incident on the CMZ) in the isotropic outflow model is 

rout ~ -77 , (16) 

R c 

where H/R ~ 0.2 — 0.3 is the geometrical aspect ratio of the 
disc (see Fig. 1 in Morris & Serabyn 1996). 

Comparing the two for the fiducial parameters accepted 
above, we have 



5 ACCRETION VERSUS STAR FORMATION 

One of the major puzzles in how SMBHs gain mass is the 
role of accretion disc self-gravity in the outer £ 0.03 — 0.1 pc 
regions. Here we expect the disc to be very cold and massive, 
provoking the conversion of gas into stars (e.g., Paczynski 
1978, Kolykhalov & Sunyaev 1980, Goodman 2003). Hy- 
drodynamical simulations of planar accretion discs in the 
regime appropriate for the Sgr A* star formation episode 
lead to rapid gas depletion through star formation (Nayak- 
shin et al. 2007), leaving almost no fuel for the SMBH. The 
Fermi-LAT observations suggest that at least in the last star 
formation episode near Sgr A* it managed to gain roughly 
the same amount of gas as was used to make stars, which 
appears to be the first observational evidence of this kind. 

The interpretation of the physical significance of an am- 
icable ~ 50% split of the gas consumption between the stel- 
lar disc and Sgr A* is model dependent. On the one hand, 
it may be an evidence that feedback from stars inside the 
accretion disc is able to stave off star formation for long 
enough to channel a sufficient amount of fuel to the SMBH 
(Thompson et al 2005). Alternatively, the dynamically hot 
structure of the young stars in the central parsec of the 
Milky Way is best explained by a non-planar gas deposi- 
tion event resulting from, e.g., collisions of two massive gas 
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clouds (Hobbs & Nayakshin 2009), collision between a cloud 
and the circumnuclear disc, or capture of a large, turbulent 
giant molecular cloud (Wardle & Yusef-Zadeh 2008, Bon- 
nell & Rice 2008). If this is so, a fraction of the gas may 
have small enough angular momentum to orbit within the 
innermost non self-gravitating disc region, avoiding the star 
formation catastrophe entirely (King & Pringle 2007, Nayak- 
shin & King 2007, Hobbs et al. 2010). In particular, in all 
of the cases simulated by Hobbs & Nayakshin (2009), the 
mass of the gas captured inside their inner boundary (a non 
self-gravitating part of the disc) was comparable with the 
mass required to fuel the Fermi-LAT lobes. 



6 DISCUSSION 

We have shown that the shape and energy content of the 
gamma-ray lobes observed by Fermi are consistent with the 
effects of a near-isotropic outflow driven by the Milky Way's 
last Eddington outburst, which we hypothesize to have been 
coincident with the well-known star formation event in the 
central 0.5 pc of the Galaxy about 6 Myr ago. The shape fol- 
lows because such outflows cannot penetrate the galaxy disc, 
in agreement with the recent conclusion that SMBHs do not 
correlate observationally with galaxy discs (Kormendy, Ben- 
der & Cornell 2011). The accreted mass driving this event is 
> 2 x 1O 3 M0, comparable with the mass of the ring of stars 
born in the event. This suggests that a significant fraction 
of gas captured into orbit by the SMBH has low enough an- 
gular momentum to accrete rather than being turned into 
stars. Because the outflow is Rayleigh-Taylor unstable, the 
gamma-ray producing gas is homogeneously mixed though 
the lobes. 

Our model is similar to that of Crocker & Aharonian 
(2011) in terms of assuming the CR protons (rather than 
electrons) produce the observed gamma-ray emission. We 
differ in the source and the age of the lobes. It is worth 
noting that if the bubbles are indeed very old quasi-static 
features as suggested by Crocker & Aharonian (2011), then 
our estimates on the recent Sgr A* activity should be taken 
as a non-trivial upper limit for that activity. 

The picture of the recent Sgr A* outburst developed 
here fits with a general view that all galaxies are active and 
reach the Eddington limit from time to time. This latter 
property is seen most obviously among galaxies classified 
on other grounds as active, particularly in the very high 
fraction of local AGN with high-speed outflows (Tombesi 
et al. 2010a, b). It seems likely that other 'normal' galaxies 
may produce similar phenomena. However the comparative 
dimness (L 7 ~ 10 37 ergs -1 ) of the lobes may make these 
difficult to detect. A prediction of the present paper is that 
the outermost parts of the lobes should still be expanding 
at ~ 1000 kms -1 . 
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